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ABSTRACT 

Context. HD 181327 is a young Main Sequence F5/F6 V star belonging to the f? Piet, oris moving group (age ~ 12 Myr). It harbors an 
optically thin belt of circumstellar material at ~90 AU, presumed to result from collisions in a population of unseen planetesimals. 
Aims. We aim to study the dust properties in the belt in great details, and to constrain the gas-to-dust; ratio. 

Methods. We obtained far-IR photometric observations of HD 181327 with the PACS instrument onboard the Herschel Space 
Observatory*, complemented by new 3.2 mm observations carried with the ATCA** array. The geometry of the belt is constrained 
with newly reduced HST/NICMOS scattered light images that break the degeneracy between the disk geometry and the dust 
properties. We then use the radiative transfer code GRaTer to compute a large grid of dust models, and we apply a Bayesian 
inference method to identify the grain models that best reproduce the SED. We attempt to detect the oxygen and ionized carbon 
fine-structure lines with Hersehel/PACS spectroscopy, providing observables to our photochemical code ProDiMo. 

Results. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected with Herschel/PACS 
completing nicely the SED in the far-infrared. The disk is marginally resolved with both PACS and ATCA. A medium integration 
of the gas spectral lines only provides upper limits on the [OI] and [CII] line fluxes. We show that the HD 181327 dust disk consists 
of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an important, layer of ice, for a 
total dust mass of ~ 0.05A/®. We discuss evidences that the grains consists of fluffy aggregates.The upper limits on the gas atomic 
lines do not, provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than ~ 17 Ala. 
Conclusions. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that 
it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive 
Edgeworth-Kuiper Belt, that may be a source for the future delivery of waf er and volatiles onto forming terrestrial planets. 

Key words, stars: individual: HD 181327 circumstellar matter infrared: planetary systems radiative transfer 


1. Introduction 

The evolution of planets in the Solar System is intimately 
connected to the existence of a reservoir of planetesimals 
in its outer regions. The depletion of the Kuiper Belt and 
the main asteroid belt that occurred during the Late Heavy 
Bombardment not only delivered large amounts of water, 
volatiles and carbonaceous material onto the inner plan- 
ets, but it also reduced the later rate of catastrophic im- 
pacts onto the Earth, opening the path to the emergence 
of life (Morbidelli et, al. 2000; Greaves & Wyatt 2010; 
Hartogh et al. 2011). Although detecting Kuiper Belt-like 
objects around nearby planetary systems will remain an 
unachievable goal for the foreseeable future, their colli- 
sional erosion produces circumstellar disks of dust respon- 
sible for a characteristic excess emission detectable at in- 
frared and sub-millimeter wavelengths. Recent, studies us- 
ing the Herschel Space Observatory reveal that as much 
as ~30% of nearby Main-Sequence F, G. K stars are sur- 


Send offprint requests to: jeremy.kbretonQobs.ujf-grenoble.fr 
* Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consor- 
tia and with important participation from NASA. 

** The Australian Telescope Compact Array is operated by 
the Australian Telescope National Facility (ATNF) managed by 
CSIR.O. 


rounded by cold debris rings analogous to the Kuiper Belt 
(latest, Hersehel/DUNES open time key program results, 
Eiroa et al. in prep). However, the cold temperature and 
the large grain sizes make it impossible to detect solid-state 
features in their spectrum, thus preventing unambiguous 
identification of the dust composition. We are thus left with 
model-dependent methods to constrain the dust composi- 
tion from color and polarimetric measurements, and from 
the Spectral Energy Distribution (SED). 

Another uncertainty in planet formation theories lies 
in the lack of detailed observational constraints for the 
mutual dust and gas dissipation timescale in circumstel- 
lar disks, or more precisely, on how the gas-to-dust, mass 
ratio evolves with time and location in disks at, the very 
early stages of planet formation. The study of stellar clus- 
ters of different ages shows that, statistically, a population 
of young stars loses its massive dust disks in only a few 
million years on route to the Main Sequence (Carpenter 
et al. 2005), although, for individual objects, the transi- 
tion from an optically thick to an optically thin disk is 
expected to be even quicker (a million years or less, Cieza 
et al. 2007; Currie & Sieilia-Aguilar 2011). But not much 
is known about the characteristic evolutionary timescales 
for the gas in these disks. Observing young stars at dif- 
ferent evolutionary stages can help elucidate the extent to 
which gas and dust simultaneously dissipate in disks. This 
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requires observations of both the dust and gas components 
and, at the same time, models able to deal with both the 
continuum and line emission. 

As part of the GASPS (Gas in Protoplanetary 
Systems) Herschel Open Time Key Programme, we ob- 
served IID 181327, a young (~12 Mvr) F5.5V star located 
at 51.8 pc (Holmberg et al. 2009, Tab. 1) and belonging to 
the 8 Pictoris moving group (Zuckerman & Song 2004; 
Mamajek et al. 2004). Its Spectral Energy Distribution 
identifies HD 181327 as a main sequence debris disk candi- 
date with a fractional infrared luminosity Tir/jL* = 0.2% 
(Mannings & Barlow 1998). Schneider et al. (2006), from 
NICMOS coronagraphic observations in scattered light, dis- 
covered a ring-like disk of circumstellar dust located at ~ 
89 AU from its star. The ring is inclined by 31.7+1.6° from 
face-on and shows an apparent brightness assymetry with 
respect to the minor axis that is well explained by a strong 
directionally preferential scattering. No photocentric offset 
is seen in the ring relative to the position of the central star 
and the disk appears axisymmetrical. Recent reprocessing 
of these data (Schneider et al.. in prep.) confirms the an- 
nulus shape of the disk, and suggests a narrow radial dis- 
tribution of the dust in the system (FWHM = 24.5 AU). A 
low surface brightness diffuse halo is seen in the NICMOS 
image and in complementary HST/ACS images at a dis- 
tance of ~ 4" < r < 9" from the central star which may 
correspond to a population of barely bound grains origi- 
nating in the main ring. Mid-infrared images at 18.3 pm 
by Chen et al. (2008) suggest asymmetries in the density 
profile of the main ring with respect to the apparent mi- 
nor axis. These 18.3pm observations also show that little 
to no emission (< 15%) comes from inside the ring re- 
solved in scattered light. Therefore, the density profile of 
the HD 181327 dust disk should be well constrained by the 
scattered light and thermal emission images, allowing a de- 
tailed study of the dust properties. Interestingly, from IRS 
and MIPS-SED data, Chen et al. (2008) also concluded on 
the presence of ice in its crystalline form in the debris disk. 

In the present paper we perform spectral energy 
distribution (hereafter SED) and line fitting of new 
Herschel/PACS observations of HD 181327 (presented in 
Sec. 2), supplemented by data from the literature. We also 
present new ATCA imaging of the disk (Sec. 2). In Sec. 3, 
we derive surface density profiles from the newly reduced 
HST /NICMOS data, which well constrain the ring geome- 
try. This alleviates the degeneracy between disk structure 
and dust properties in debris disk modeling, and allows us 
to perform a detailed study of the composition and size 
distribution of the dust grains. We describe our dust model 
and present the results of the SED fitting in Sec. 3 and 
4. We assess the gas content in this young debris disk in 
Sec. 5. Finally we discuss the results in Sec. 6, 7 and 8 and 
present our conclusions in Sec. 9. 

2. Observations and Data Reduction 

2.1. Herschel/PACS continuum 

We obtained new far-IR photometry and spectroscopy of 
HD 181327 using the PACS instrument on the Herschel 
Space Observatory (Pilbratt et al. 2010; Poglitsch et al. 
2010). We begin by describing the photometry observa- 
tions and data reduction. The observations were carried 
out in two modes: 1) 159 sec of point-source chop-nod mode 


Table 1. Star and disk properties 


Parameter 

Value 

Reference 

Spectral Type 

F5/F6V 

Nordstrom et al. (2004) 

Magnitude V’ 

7.0 

Bessell (2000) 

Gravity (log(g)) 

4.510 

Chen et al. (2006) 

Distance (d*) 

51.8 pc 

Holmberg et al, (2009) 

A „ 

12 Mvr 

Zuckerman & Song (2004) 

Age 


Mamajek et al. (2004) 

Luminosity (L*) 

3.33 1$ 

this study 

Mass (Mj) 

1.36 M© 

this study 

Disk PA' b> 

1074:2° 

Schneider et al. (2006) 

Disk inclination 

31.7±1.6° 

Schneider et al. (2006) 

Tir/T* 

~ 2 x 10~ 2 3 

this study 

<a) From integrating the synthetk 

; stellar spectrum, ^'Position 


Angle, East of North. 

imaging at 70 and 160 pm (obsid: 1342183658) and 2) two 
276 sec scan map imaging observations at 100 and 160 pm 
(obsids: 1342209057-8). The scan maps were executed with 
the medium scan speed (20"/s), at scan angles of 70° and 
110°. The two cross scan maps at different scan angles were 
combined to increase the sigrial-to-noise by cutting down on 
streaking in the image background along the scan direction. 

We reduced the photometry with the Herschel 
Interactive Processing Environment (HIPE, Ott 2010) ver- 
sion 4.2. The uncertainties in the absolute flux calibration 
for this version are 10% at 70 and 100 pm and 20% at 
160 pm (see the PACS Scan Map AOT release note 1 ). For 
the 70’’ and 160 pm chop-node images, we performed aper- 
ture 1 photometry using apertures of 7" and 11", respectively. 
For the 100 and 160 pm scan map observations, we used 
apertures of radius 12" and 24", respectively. An annulus 
for sky background estimation was placed from 10" - 20" 
beyond the aperture in each case, and an aperture correc- 
tion was applied to the measured fluxes. Both estimations 
at 160 pm are fully compatible within the error bars; the 
ICO pm scan-map flux is used in the rest of the study. Since 
the source is marginally resolved at 70 and 100 pm (see 
Section 2.1), we calculated our own aperture corrections 
following the procedure described in Poglitsch et al. (2010). 
The photometry results are presented in Table 2. 

Table 2. Photometry results for HD 181327 from chop-nod a> 
and scan map - 2 ' observations. 


Wavelength 

(pm) 

F f 

cont 

tty) 

Absolute Flux Cal. 
U ncertainty 

70 (I) 

1.827 ±0.0069 

10% 

100 f2) 

1.337 ±0.0082 

10% 

160 (1) 

0.767 ± 0.015 

20% 

160 (2> 

0.772 ±0.011 

20% 


’The +1<T error bars on the flux measurements include 
statistical uncertainty only. 


HD 181327 is more extended than the Herschel point- 
spread function (PSF) in both the 70 pm and 100 pm im- 
ages. The 70 pm image of HD 181327 appears in Figure 1, 
along with the 70 pm image of a PSF reference star ob- 
served using the same chop-nod mode (HD 148387). The 

1 available at http://herschel.esac.esa.int/ 

AQTsReleaseStatus . shtml 
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top panel of Figure 2 shows the azimuthally-averaged ra- 
dial brightness profiles for HD 181327 and HD 148387. For 
the 100 //in scan map observation of HD 181327, we used 
a different PSF reference observation, i.e. a scan map of 
a Boo. The 100 /on a Boo radial profile appears in Figure 2, 
bottom panel. At both wavelengths, the radial profile of 
HD 181327 is clearly extended beyond that of the PSF ref- 
erence. The FWHM of the HD 18i327 70 //in profile is 6744, 
compared to a FWHM of 5"28 for HD 148387. At 100 //in, 
the FWHM for FID 181327 is 7758, compared to 6794 for 
a Boo. 


2.2. ATCA 3.2mm observations 

The Australia Telescope Compact Array (ATCA) was used 
to observe HD 181327 at 3.2 mm in May and September 
2010. The observations were carried out in continuum dou- 
ble sideband mode, with 2046 x 1 MHz channels per side- 
band with an effective total bandwidth of just over 2 GHz. 
The ATCA has five antenna with 3-mm receivers, and our 
observations were carried out with the compact H75 ar- 
ray configuration with baselines of 31 98 m. The quasars 
PKS 1815-553 and PKS 1253-055 were used for the gain and 
bandpass calibrations respectively, and Uranus was used for 
the absolute flux calibration. The ATCA calibration uncer- 
tainty is estimated to be about. 20% in the 3-mm band. 
The data were calibrated and reduced with the MIRIAD 
package (Sault et al. 1995). The frequency pair for the ob- 
servations were centered at 93 and 95 GHz, and during 
the reduction process the data was averaged over 64 chan- 
nels to improve the signal-to-noise, shifting the frequency 
pairs to 92.968 and 94.969 GHz. The final CLEANed map 
combined the two sidebands, resulting in an image at 
93.969 GHz (3192.5 //m). Our total on-source integration 
time was 770 minutes. We determined an integrated flux of 
7.2 x 10 -4 Jy, with an RMS of o — 1.1 x 10~ 4 Jy/beam. 
Including the 20% calibration uncertainty, this yields an 
uncertainty of 2.5 x 10~ 4 Jy. The disk is marginally resolved. 
Assuming a Gaussian profile, we find a PA of ~ 102°, con- 
sistent with that seen on the NICMOS image (107 ± 2°. 
Tab. 1), and along which the FWHM is ~ 7.8" (beam: 
3.5" x 2.8", PA = 85.4°). The center of the Gaussian fit, 
is also shifted by 0.6" toward the West with respect to the 
star position. Structures consistent with a belt at ~100 AU 
arc detected along the semi-minor axis with a peak flux of 



Fig. 1. PACS 70 fim imaging of (left) HD 181327 and (right) 
a PSF reference star, HD 148387. Both images were acquired in 
chop-nod mode. The field of view of each image is 25" x 25". 
Brightness contours at 5, 10. 25, 50 and 75 times the background 
noise in the HD 181327 image are overlaid on both images. 



Fig. 2. Radial profiles of HD 181327 at 70 //in (top), compared 
to HD 148387 as a PSF reference, and at 100 // m (bottom), 
compared to a Boo as a PSF reference. Both plots show the 
azimuthally averaged, normalized radial brightness profile vs. 
radius from the. star. Each error bar is the standard deviation 
of the brightness values at that radius. The FWHM values for 
the profiles are given on the plot. 


4.01 x 10 4 Jy/beam oil the NW side. A possible brightness 
asymmetry with respect to the SE side is revealed (Fig. 3). 

2.3. Herschel/PACS Spectroscopy 

Our HD 181327 spectroscopy was obtained with the PACS 
instrument, which is an IFU (Integral Field Unit) with a 
5x5 array of 974 x 974 spaxels. The data were taken in 
chop-nod mode using the line scan and range scan modes. 
The spectra cover 8 small wavelength ranges (~ 1.8 pm 
wide) targeting particular emission lines, including [O I] at 
63 pm and [CII] at 158 pm. The spectral resolution of the 
data varies between R « 1100 and 3400. Additional details 
on the observations appear in Mathews et al. (2010). 

The data were calibrated with HIPE v4.2. The 
wavelength-dependent flux and aperture corrections pro- 
vided by the PACS development team were applied, giv- 
ing an uncertainty in the absolute flux calibration of 30% 
(see the PACS Spectroscopy Performance and Calibration 
Document 1 ). Significant emission was detected only in the 
central spaxel. No emission lines are seen, only continuum 
flux. The central spaxel spectrum centered on the [O I] 
63 / mi line appears in Figure 4. 

The continua were binned at half the instrumental 
FWHM and fitted with 1st degree polynomials, using the 
fluxes and errors produced by our HIPE reduction script . 
The HIPE errors appeared to be slightly too small, however. 
New flux random errors were calculated from the RMS of 
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Fig. 3. ATCA 93.969 GHz (3192.5 mm) map of HD 181327 with 
contours at -2 (dashed red line), 2 and 3 a (solid blue lines), 
where a is the image RMS given by l.lxlO -4 Jy/beam. The 
dashed blue area represents the beam. The Gaussian fit is over- 
plot with contours at 2 and 3<r (solid green lines). The red 
and green crosses denote respectively the star position and the 
Gaussian fit center. 


2.5 



x 

Cl 1.5 


1.0 


01 63 pm 
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Fig. 4. PACS spectrum of HD 181327, centered on the wave- 
length of the [O I] 63 pm emission line. The linear fit to the 
continuum flux is overplotted with the red line. The [O I] 63 pm 
emission line was not detected. 


and the gas, as well as on the dust properties. These are de- 
tailed below in Secs. 3.1 and 3.4. Additional data collected 
from the literature and stellar properties are summarized 
in Sec. 3.2. 


the data minus the continuum fits, over wavelength ranges 
spanning 6 to 14 times the instrumental FWHM, centered 
on the expected line rest wavelengths. The continuum fluxes 
at the rest wavelengths of the expected emission lines ap- 
pear in Table 3. To determine the upper limit on the in- 
tegrated line flux for each available emission line, we inte- 
grated the data over a small wavelength range, either ± 1 or 
±1.5 resolution elements, centered on the rest wavelength 
of the expected emission line. The flux errors were propa- 
gated to determine the Icr upper limit on each integrated 
line flux. The 3 a integrated line flux upper limits appear in 
Table 3 as well. 


Table 3. Spectroscopy Results for HD 181327 


Line 

Wavelength 

(fim) 

r cont 

(Jy) 

f'L 

(xlO -18 W/m 2 ) 

[O I 

63.1852 

1.731 ± 0.236 

< 9.82 

CO 36-35 

72.8429 

1.491 ±0.185 

< 18.41 

o-H 2 0 

78.7414 

1.431 ±0.187 

< 10.51 

CO 29-28 

90.1630 

1.264 ± 0.222 

< 8.90 

[O Ij 

145.5350 

0.804 ±0.107 

< 8.48 

[CII] 

157.6790 

0.767 ±0.175 

< 7.02 

0 -H 2 O 

179.5265 

0.388 ± 0.201 

< 6.87 

DCO i ‘ 22-21 

189.5700 

0.581 ± 0.378 

< 16.71 


7 The ±1 <t error bars on the flux measurements include sta- 
tistical uncertainty only. The absolute calibration uncertainty 
amounts to 30%. 1 3 a upper limit 


3. Model Setup 

The dust and gas models developed in Secs. 4 and 5 rely 
on some assumptions on the spatial distribution of the dust 


3.1. New Surface Density Profile for the HD 181327 Ring 

In this study, we wish to take advantage of the fact that 
the spatial distribution of the dust in the HD 181327 disk 
is well-constrained, in order to examine the dust properties 
in details and to get an upper limit on the gas mass from 
the Herschel non-detections of gaseous emission lines. The 
disk around HD 181327 has indeed been spatially resolved 
in scattered light with the HST/NICMOS facility at 1.1 /jm 
by Schneider et al. (2006), which alleviates many uncertain- 
ties on the disk structure. The NICMOS images revealed 
a ring- like disk with a peak density position rp = 1.71” 
(~ 88.6 AU, according to the new distance to the star, 51.8 
pc, by Holmberg et al. 2009). with a FWHM of ~ 0.538” 
(~ 28 AU). 

Recently, Schneider et al. (2011, in preparation) im- 
proved the data reduction of the original HST/NICMOS 
observations of HD 181327 thanks to a better selection of 
the reference PSF to be subtracted from the raw data. This 
results in an overall improvement of the signal-to-noise, es- 
pecially inside the ring, and leads to a slight sharpening 
of the belt after sky plane deprojection and taking the az- 
imuthal median of the disk image. The new profiles, dis- 
played in the top panel of Fig. 5, show a peak at a dis- 
tance of 88.6 AU and a surface brightness fainter by ~ 21% 
(0.76mJy.arcsec~ 2 at the peak position) compared to the 
original Schneider et al. (2006) study. 

We derive the shape of the surface density profile £(r) 
by inverting the new scattered light surface brightness pro- 
file. We employ the methodology described in Augereau & 
Beust (2006), assuming a disk inclination of i — 31.7° from 
face-on. The inversion procedure in fact yields the product 
<r sca E(r), where er sca is the mean scattering cross section of 
the grains at 1.1 fim. This is equivalent to the product of the 
vertical optical depth r±(r), multiplied by the mean albedo 
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Distance from the star in AU 

Fig. 5. New HST/NICMOS 1.1 pm surface brightness profiles 
(top. nominal and ±lcr profiles), and surface density profiles 
resulting from direct inversion of the brightness profiles (bot- 
tom), assuming two possible values for the asymmetry param- 
eter. The color-coding is the same for both panels. The striped 
areas identify regions dominated by PSF-substraction residuals 
where the, density has been extrapolated . 


u ) at that wavelength (cr &ca E(r) = wt±(t), see Eqs. 2-3 in 
Schneider et al. 2006). The solution of the inversion proce- 
dure depends on the anisotropic scattering phase function, 
and on the disk vertical profile. Nervertheless, Schneider 
et al. (2006) showed that the density profile is largely in- 
dependent of the adopted shape for the vertical dust dis- 
tribution in the range of acceptable vertical scale heights 
Hq. We here assume a Gaussian vertical profile with a scale 
height II q = 6 AU at the observed peak radius of the ring 
To = 88.6 AU (Ho/ro ~ 0.068, consistent with the theoreti- 
cal lower limit of 0.04 from Thebault 2009). The observed 
phase function of the disk at 1.1 pm is well-represented by a 
Henyey-Greenstein scattering phase function with asymme- 
try parameter |</hg|= 0.30 ±0.03. This value was discussed 
by Chen et al. (2008) in light of asymmetries identified in T- 
ReCS thermal emission images in the Q„-band. Chen et al. 
(2008) suggest that the asymmetries may have led to an 
overestimate of |.9 hg|, which may be as small as 0.16±0.04. 
We thus perform the inversion for both |c/hg |== 0.16 and 
|.9 hg|= 0.30. 

The results are displayed in the bottom panel of Fig. 5, 
where we also show the impact of the surface brightness un- 
certainties on the inferred surface density profile. We find 
an actual surface density peak at 89.5 AU from the star. 
A higher |<?hg| leads to a globally higher surface density 
profile. The surface brightness measured interior to ~ 1.5” 
(~ 78 AU) is not reliable because of PSF subtraction resid- 
uals. We lit the profiles from 79 to 83 AU obtained for both 
!<7hg | values by a very steep power-law function of the form 
~ r +10±1 . We extrapolate this profile interior to 78 AU as 
shown in Fig. 5 (striped area). This extrapolation is sup- 
ported by the mid-IR results of Chen et al. (2008), who 


show that the inner part, of the disk is largely cleared to- 
ward the center, containing little to no dust. 

Outside of the belt peak position, the improved data 
reduction results in a profile that falls off as ~ r -4.7±o.3 
from ~92 to 162 AU, a steeper fall-off than the Schneider 
et al. (2006) original profiles (~ r~ 3 ). The dust is thus 
more tightly confined to the parent-belt, than previously 
estimated, and the profile is also steeper than theoretical 
expectations. The exact value of the outer slope can only 
slightly impact the modelling results, as most of the emis- 
sion arises from the inner part of the belt where the temper- 
ature is larger, except possibly for the longer-wavelengths. 

Eventually, the differences between the several radial 
profiles in Fig. 5 essentially translate into a scaling factor 
that will impact the total disk mass estimate. We estimate 
that the uncertainty on the asymmetry parameter adds an 
additional ~4% to the uncertainty on our final mass values. 
The density profile obtained for |,9hg |=0.30 is used in the 
rest of this study. We stress that, in order to focus on the 
other disk properties, we make the assumption that this 
profile is valid for any dust, grain constrained by the SED, 
i.e. for grain sizes ranging from the micron to the millimetre 
scale. This assumption remains reasonable as long as the 
flux at any wavelength arises predominantly from regions 
close to the density peak position (see Sec. 6.2). 


3.2. Spectral Energy Distribution 

The Herschel and ATCA data presented in this paper are 
complemented by IR and sub -111111 data collected from the 
literature. These are listed in Table 4. The shortest, wave- 
length flux we use is the 12 pm IRAS photometry for which 
the emission shows 110 excess from the stellar photosphere 
level. This observation is documented in Schneider et al. 
(2006) together with IR-excess from the IRAS Faint Source 
Catalogue and obtained with the Spitzer/MIPS instrument. 
These are complemented by ground-based mid-IR measure- 
ments performed with the T-ReCS imager on Gemini-S by 

Chen et al. (2008), and 13 points extracted from the 5 

35 pm Spitzer/IRS spectrum that we reduced using the c2d 
Legacy team pipeline (Lahuis et, al. 2006). 

Additional MIPS observations were obtained in SED- 
mode by Chen et al. (2006) between 54.5 and 95.5 pm, 
well-constraining the shape of t he disk SED around its max- 
imum emission wavelength. Far-IR excess were reported by 
Mooret al. (2006) using the ISO database. The shape of the 
SED in the far-IR domain (100-200 pm) is well constrained 
by our new Herschel measurements. 

Because the uncertainties on the dust optical constant s 
at millimetre wavelengths are high, we chose not to use 
the ATCA flux to constrain the dust model (see Sec. 6.2 
for a discussion on the compatibility of the model to these 
observations). Therefore, the sub- mm regime is only rep- 
resented by a measurement with the LABOCA bolometer 
array at, the 12-m telescope APEX (Nilsson et al. 2009, 
A ~ 870 pm), which thus represents a key piece of informa- 
tion. These archival data, together with our new observa- 
tions, constitute the 49 measurements that are displayed in 
Fig. 6. 
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Table 4. HD 181327 photometric observations from the litera- 
ture. The Herschel continuum values used for the dust modeling 
are in Tables 2 and 3. 


Wavelength 

(pm) 

Flux 

(Jy) 

Uncertaintv 

(Jy) 

Instrument 
k Reference 

12 

0.114 

0.011 

IRAS 

[1J 

15.0 

0.094 

0.004 

Spitzer/IRS 

[2] 

16.3 

0.097 

0.005 

Spitzer/IRS 

[2| 

17.4 

0.107 

0.005 

Spitzer/IRS 

|2| 

18.3 

0.114 

0.003 

T-ReCS 

13) 

18.7 

0.125 

0.007 

Spitzer/IRS 

12] 

19.8 

0.144 

0.009 

Spitzer/IRS 

|2| 

20.4 

0.149 

0.008 

Spitzer/IRS 

|2| 

21.1 

0.157 

0.020 

Spitzer/IRS 

121 

23.1 

0.224 

0.007 

Spitzer/IRS 

|2| 

23.7 

0.223 

0.009 

Spitzer/MIPS 

HI 

25 

0.306 

0.024 

IRAS 

)1| 

25.7 

0.303 

0.014 

Spitzer/IRS 

[21 

28.0 

0.384 

0.014 

Spitzer/IRS 

2| 

30.6 

0.488 

0.026 

Spitzer/IRS 

|2| 

33.1 

0.623 

0.025 

Spitzer/IRS 

[2] 

35.5 

0.718 

0.034 

Spitzer/IRS 

|2| 

54.5 

1.88 

0.09 

MIPS-SED 

[4) 

58.0 

1.83 

0.07 

MIPS-SED 

[4] 

60 

1.73 

0.17 

ISO 

[3] 

61.4 

1.98 

0.08 

MIPS-SED 

[4] 

64.9 

2.01 

0.08 

MIPS-SED 

r 4] 

68.4 

1.98 

0.08 

MIPS-SED 

[4] 

71.4 

1.73 

0.12 

Spitzer/MIPS 

HI 

71.6 

1.99 

0.09 

MIPS-SED 

4J 

75.1 

1.93 

0.10 

MIPS-SED 

[4] 

78.4 

1.81 

0.10 

MIPS-SED 

'4] 

81.8 

1.72 

0.11 

MIPS-SED 

[4] 

85.3 

1.62 

0.11 

MIPS-SED 

[4] 

88.8 

1.47 

0.11 

MIPS-SED 

4] 

90 

1.41 

0.14 

ISO 

[5] 

92.0 

1.46 

0.13 

MIPS-SED 

A] 

95.5 

1.37 

0.17 

MIPS-SED 

A] 

100 

1.7 

0.2 

IRAS 

[1] 

155.9 

0.77 

0.09 

Spitzer/MIPS 

4] 

170 

0.736 

0.192 

ISO 

.5) 

870 

0.052 

0.006 

LABOCA 

[6] 

3190 

7.2xl0“ 4 

2.5xl0 -4 

ATCA 

A] 


Notes and References - All uncertainties are 1<t and include 
absolute calibration error. [1] Schneider et al. (2006), [2] This 
study, data reduced with the Spitzer/e2d team pipeline (Lahuis 
et ai. 2006), [3] Gemini, Chen et al. (2008), [4] Spitzer, Chen 
et al. (2006), [5] Moor et al. (2006). [6] APEX. Nilsson et al. 
(2009), [7] This study, Sec. 2.2. 


3.3. Star Properties 

We derive excesses from a comparison of the measure- 
ments with a synthetic stellar spectrum extracted from 
the “NextGen” model atmosphere grid (Hauschildt et al. 
1999) and scaled to match the visible magnitude V = 7.0 
mag (Bessell 2000), assuming Iog(<?)=4.5, T e g — 6600K, 
which leads to a stellar luminosity of 3.33 Lq and a mass 
of 1.36 Mg. This stellar model shows good agreement with 
the Spitzer/IRS spectrum at wavelengths shorter than 10 
15/iin. longward of which the excess starts dominating 
over the stellar photospheric emission. The star properties 
adopted in this study are summarized in Tablet. 



Fig. 6. Observed Spectral Energy Distribution of HD 181327 
(see Tables 2, 3 and 4) 


3.4. Dust Properties: Size Distribution and Composition 

For both the dust and gas modeling, we make the assump- 
tion that the grain population has a differential size dis- 
tribution valid anywhere in the ring. This is a central hy- 
pothesis that necessarily limits the uniqueness of the solu- 
tion. We discuss its impact in Sec. 6. The differential size 
distribution follows a classical power law dn(a) cx a K da, 
from a minimum size a m i„ to a maximum size a max . The 
minimum grain size a m ;„ and the power law exponent k 
are free parameters that are constrained by the SED fit- 
ting (Sec. 4). The maximum grain size a max is a fixed pa- 
rameter, sufficiently large such that it does not affect the 
light emission properties at the wavelengths considered (we 
chose a max ~ 8 mm). For standard k values (close to -3.5 or 
smaller), the mass of solid material in a debris disk can be 
’’arbitrarily” large depending on the assumed higher cut- 
off value of the size distribution (that can extend up to 
planetesimal sizes). In this article, we define the dust mass 
A/dust as the total mass of material contained in grains 
smaller than 1 mm. This definition has been followed within 
the GASPS consortium because it allows direct comparison 
between different studies independently from the model as- 
sumptions. 

We elaborate a complex grain model that bears some 
similarities to those succesfully employed by several authors 
for fitting SEDs of young debris disks (Li &t Greenberg 
1998; Augereau et al. 1999; Li & Limine 2003a, b), and 
for fitting scattered light disk colors (Kohler et al. 2008). 
We consider icy porous aggregates, with a silicate core 
and an organic refractory mantle (or simplified versions 
of this model; see Sec. 4), to mimic fluffy particles made 
of essentially unalterred interstellar material (the so-called 
“cold-coagulation” dust model in Li & Limine 2003a, b). 
Assuming cosmic abundances, the volume of the carbona- 
ceous mantle Vc is estimated to be ~ 2.12 times that of the 
silicate core 1% in the interstellar medium, a value similar 
to the one derived from the mass spectra of comet Halley 
in the Solar System (Vfc/V si ~ 2, Greenberg 1998). 

To calculate the optical properties, we use a num- 
ber of representative materials for which optical indexes 
are tabulated. These are either derived from laboratory 
measurements or from observational and analytical argu- 



J. Lebreton et al.: An icy Kuiper-Belt around HD 181327 


7 


Table 5. Parameter space explored with GRaTer , corresponding to a total of about 1,150,000 models. The volume fractions t'si. 
vc. Vice and the porosity V are defined in Sec. 3.4 (Eq. 1) 



Description 

Parameter 

Explored range 

Values 

Distribution 


minimum grain size 

Oinin [pm] 

0.01... 100 

77 

log 


size power law index 

K 

-6.0. . . -2.5 

50 

linear 


maximum grain size 

Umax [mm] 

7.8 

fixed 



Disk mass up to 1 mrn 

.Ui„M [AL®] 

> 0 

fitted 


1 AND 2-MATERIAL MODELS 

Sec. 4-1 ■ 1 

Si+ Amorphous C 

ACAR volume fraction t>c 

0.00. . . 1.00 

21 

linear 

Sec. 4-1.2 

Si+amorphous ice 

ice volume fraction thee 

0.00 ...0.90 

19 

linear 

Sec. 4-1-2 

Si+crystaliine ice 

ice volume fraction v tC e 

0.00 . . . 0.90 

19 

linear 

Sec. 4-1-3 

Porous Si 

Porosity V 

0.00 . . . 0.99375 (a) 

24 

linear 

3 AND 4-material models 


Si 

astrosi volume fraction vs. 

l/3x (l-l>jce ) 



Sec. 4-2 

+amorphous C 

am. carbon volume fraction vc 

2/3 X (1-Uj ce ) 



-(-amorphous ice 

ice volume fraction vi ce 

0.00. . . 0.90 (b) 

12 

linear 


+■ vacuum 

Porosity V 

0.00. . . 0.95 <c) 

18 

linear 


Notes (a) 19 linearly distributed points from 0.00 to 0.95 + additional high porosity points : 0.96875, 0.975, 0.9875, 0.99375. 
99 10 linearly distributed points from 0.00 to 0.90 + additional points at 0.65 and 0.75. (c) 10 linearly distributed points from 0.00 
to 0.90 + additional points at 0.45, 0.55, 0.625, 0.65, 0.675, 0.75, 0.85, 0.95. 


1 000.0 E 


100.01 


S' 

5 loot 


1-Ot 


o.i L 
0.- 


Compact Si (p“3.50g/cm\ ^ 2 =10.3) 
Si+amorph. ice (p=2.10g/cm 3 , % 2 =8 30) 
Si+cryst. Ice (p=2.27g/cm 3 , x 2 =10.4) 
Amorph. carbon (p=2.20g/cm 3 , x 2= 7.75) 
Highly porous Si (p=0.09g/cm 3 . x 2 =1 1 .9) 
Si+C+ice+Poros. (p=0,16g/cm 3 , y 3 =1.64) 


'■-•IV 

•v; 


1.0 


10.0 
X [pm] 


\ 


100.0 


1000.0 


Fig. 7. Expect gas line fluxes as a function of gas-to-solid mass 
ratio for either a low-PAH abundance (thin lines) or a high-PAH 
abundance (thick lines). Horizontal lines denotes the Herschel 3a 
upper-limits. The gas to solid mass ratios assume a solid mass of 
0.164 M®. The PACS 3-cr upper limits do not include the 30% 
flux calibration uncertainty. 


been shown to be a good approximation for the absorption 
and scattering efficiencies, provided that grain shapes are 
close to ellipsoids such that the porosity can effectively be 
seen as the volume fraction of vacuum (Sheri et al. 2009, 
and references therein). In fact, Shen et al. (2009) show 
that only the polarization and the phase functions are sig- 
nificantly sensitive to grain shapes. We define the composi- 
tion of a grain through the volume fractions of its different 
components : 


Vsi 

v si = 77 — , v c 

*rn at 


J£c 

V m 


“» ^ice 


Vic 

Vrnat 


Kac 

Vmat + Krac 


(1) 


where Vsi, Vc, V ce and Kac are the total volume of astro- 
nomical silicates, carbonaceous material, ice and vacuum 
in the grain, respectively, and V mat = Vsi + Vc + Vice is 
(lie total volume of solid material. V denotes the poros- 
ity. Should one of the materials within the grains reach its 
sublimation temperature, it is replaced by vacuum and the 
optical properties are adjusted accordingly in the model. In 
most cases however, the disk is too cold for any material 
to sublimate. Fig. 3.4 shows the mass opacity (absorption 
cross section per unit mass) calculated using Eq. 4 of Draine 
(2006) for several prototypical dust compositions, with the 
size distributions described in Tab. 6. 


ments. The materials are: (1) the astronomical silicates (as- 
troSi) from Draine (2003), a theoretical material consist- 
ing of amorphous silicates polluted by carbonaceous ma- 
terials and other metals (p = 3.5 g.cm -3 ); (2) the ACAR 
sample of amorphous carbon from Zubko et al. (1996) 
as organic refractories (p = 1.95 g.cm” 3 ); (3) the H ■_>()- 
dominated ice from Li & Greenberg (1998) (amorphous 
ice, p = 1.2 g.cm” 3 ); and (4) the crystalline water ice from 
Warren (1984) ( p = 0.92 g.cm” 3 ). 

The optical indexes of an aggregate are obtained using a 
standard effective medium theory (EMT), the Bruggeman 
mixing rule (Bohren & Huffman 1983), and following the 
methodology of Li & Greenberg (1998) and Augereau et al. 
(1999). The optical properties are computed with Mie the- 
ory valid for hard spheres. The EMT-Mie combination has 


4. Modeling the Dust Continuum Emission 

We use the GRaTer code (Augereau et al. 1999) to com- 
pute a grid of models that aims to reproduce the HD 181327 
Spectral Energy Distribution. GRaTer is an SED, image 
and interferometric data fitter, specially designed to effi- 
ciently model optically thin disks with parametric grain 
size and radial distributions, or distributions from dynam- 
ical simulations. The code accounts for both the scattered 
light emission and the continuum emission of dust grains 
in thermal equilibrium with a star. 

In t he present case, the spatial dust distribution is con- 
strained by the resolved HST images (Sec. 3.1, Fig. 5). The 
surface density profile derived from inverting the observed 
surface brightness was imposed in the model. Therefore, our 
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Fig. 9. Bayesian probabilities for the models in Table 6: minimum grain size, size distribution power law index, and total dust 
mass. 




Wavelength fprnj 

Fig. 8. Best fit models (Table 6) for several 1- or 2-material 
mixtures (top) and 3- to 4-material mixtures (bottom). Red 
crosses: photometric data (Fig. 6). Solid grey line: synthetic stel- 
lar spectrum. 

study of the disk continuum emission focuses on the proper- 
ties of the dust particles, namely their size distribution and 
composition. To that end, we computed large grids of mod- 
els reflecting various possible dust properties, starting from 
the simplest composition (astronomical silicates) and pro- 
gressively increasing the model complexity until obtaining 
a good match to the SED. For each model, the dust mass 
Afdust and the quality of the fit are obtained from a lin- 
ear least-square fit. between the observed and the synthetic 
SEDs. The x~ listed in this section are obtained using the 


48 data points presented in Tables 2, 3 and 4, and with the 
numbers of degrees of freedom listed in Table 6. We sub- 
sequently use a statistical (Bayesian) inference method to 
derive the probability distribution of each parameter (see 
Appendix A). 

In total we consider ~ 300 compositions, with 77 x 50 
grain size distributions for each, leading to ~ 1,150,000 
models. All the parameters of the grain models are sum- 
marized in Table 5. 


4.1. First Steps Toward a Complex Dust Model 

We start by fitting the HD 181327 SED with a single- or 
two-material grain model, which show the need for a more 
complex dust model. We consider several mixtures of the 
materials presented in Sec. 3.4 to assess the broad proper- 
ties of the grains. The model then includes 3 or 4 free pa- 
rameters: (1) the minimum grain size « min ; (2) the slope of 
the size distribution re; (3) the total mass of the disk A/dusti 
and in the case of 2-component grains: (4) the volume frac- 
tion of either carbon («c) or ice (i> ice ), or the porosity (V). 
The detailed results of the investigation are summarized in 
Table 6, the best-fitting SEDs are displayed in Fig. 8 (top 
panel) and probability curves are shown in Fig. 4. 

4.1.1. Silicates or Organic Refractories? 

The most basic grain model to consider is the case of com- 
pact spheres made of pure astronomical silicates (astroSi), 
a configuration that is often used by default in many model- 
ing approaches. We let a m - m , re and Md us t vary and we com- 
pute a grid of x 2 on which we apply a Bayesian analysis. 
We find that pure astronomical silicate grains would have 
a m j D — 0.99 ± 0.13 pm and re = —3.43 ± 0.04, leading to a 
total dust mass M c i ust = 0.35 ± 0.03A/i for grains smaller 
than 1mm (Tab. 6). These values are consistent with the 
ones typically found for debris disks, but they actually lead 
to very poor fits, as illustrated in Figure 8 (top panel) where 
the model with the smallest \ 2 (\‘ = 10.3) is displayed. 
Although the model is fairly close to the data over the 
full spectral range and agrees well with the sub-nun mea- 
surement, it fails to reproduce the shape of the emission 
peak at ~ 50 — 100 pm and the mid-IR fluxes, shifting 
the dominant part of the disk to colder temperatures. It 
overestimates by 30-40% the excess at far-IR. wavelengths 
(A ~ 100 — 100 /an), in particular the Herschel/PACS mea- 
surements. 
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Table 6. Best-fitting (Best, i.e. lowest \r) and most probable (Proha, i.e. highest probability) parameters assuming different grain 
compositions (see text for the definition of the parameters). The most probable parameters are given with 3-<r uncertainties. The 
volume fractions marked with the f symbol are imposed knowing the other free volume fractions. 


Composition 

vsi 

VC 

Vice 

V 

a m in(juni) 

K 

Afdust(M®) 

xi 

1-MATERIAL (dof = 45) 

Proba : 
Sl Best : 

1.00 




0.99T0.13 

-3.43±0.04 

0.35±0.03 


1 .00 




1.00 

-3.42 

0.35 

10.3 

„ Proba : 


1.00 



1.62±0.03 

-3.65i0.07 

0.15±0.03 


Best : 


1.00 



1.62 

-3.64 

0.15 

7.75 


2-material (dof = 44) 


Si+C w 

Proba 

0.007f 0.99±0.04 


1 ,62±0.03 

-3.65±0.07 

0.15±0.03 


Si+am.ice 

Proba 

Best 

0.59t 

0.60 

0.41 ±0.05 
0.40 

0.15±0.13 

0.16 

-3.36i0.07 

-3.35 

0.26±0.05 

0.26 

8.90 

Si+cr.ice 

Proba 

Best 

0.67f 

0.65 

0.33+0.12 

0.35 

0.29±0.24 

0.26 

-3.35±0.02 

-3.35 

0.27±0.04 

0.27 

10.4 

Si-KP 

Proba 

Best 

0.94f 

0.025 


0.057±0.062 0.98±0.20 

0.975 0.010 

-3.42±0.05 

-3.28 

0.35±0.04 

0.039 

11.9 


3-matkrTal (dof = 44) 


SiiCiP 

Proba : 
Best, : 

0.33f 

0.33f 

0.67f 

0.67f 


0.95±6xl0~ f> 

0.95 

9.45±3.20 

10.0 

-3.18±0.09 

-3.14 

0.042±0.001 

0.042 

6.05 

SiiCiice 

Proba : 
Best : 

-1 1— 

r— t- 
o o 
o © 

0.1 3f 
0.13f 

0.700±0.002 

0.80 


0.78±0.04 

0.55 

-3.50±0.08 

-3.42 

O.lliO.Ol 

0.11 

4.48 

4-material (dof = 43) 

SiiCiiceiP 

Proba : 
Best : 

O.llf 

0.12f 

0.22f 

0.23f 

0.67±0.07 

0.65 

0.63±0.21 

0.65 

0.81±0.31 

0.89 

-3.41 ±0.09 
-3.42 

0.051±0.016 

0.048 

1.64 

Notes “dof” means degrees of freedom. (1) The best fit for the Si+C mixture 

parameters for that solution are given in the 1 -material part of the table. 

is obtained for silicate-free 

grains. The best-fitting 


We take a further step and consider mixtures of amor- 
phous silicates and carbonaceous materials. We ran a 
Bayesian analysis and find that the actual best composi- 
tion corresponds to pure carbon. Qualitatively, the effect 
on the SED is not significant (see Fig. 8, top panel), even 
if it leads to a better minimum chi-square (y 2 — 7.75). 
This model gets closer to the measured fluxes at the peak 
emission (A ~ 50 — 100 ftin), but does not provide an over- 
all good fit. The addition of carbonaceous material slightly 
shifts the best fit parameters of the grain size distribution 
(Tab. 6), resulting in slightly larger a m ; n values and slightly 
steeper size distributions 

4.1.2. Icy Grains? 

Chen et al. (2008) suggested that HD 181327 may contain 
a significant fraction of water in the form of crystalline ice, 
producing a broad peak in the combined Spitzer/IRS (5- 
35 pm) and MIPS-SED (60-77 pm) spectrum. To assess the 
presence of water ice in the debris disk of HD 181327, we 
model the disk with grains made of a mixture of astronom- 
ical silicates and water ice, either in the form of amorphous 
ice or crystalline ice. 

The results of the statistical analysis are documented 
in Table 6. The quality of the fit is only slightly improved 
compared to the pure silicate model: the best models have 
Xr = 8.90 for amorphous ice, yf = 10.4 for crystalline ice. 
In fact, the addition of water ice, crystalline or amorphous, 
does not modify the general shape of the synthetic SED (see 
Fig, 8, top panel). The smallest y 2 calculated with amor- 
phous ice is slightly smaller than the one of crystalline ice, 
but in both cases the values are too high to accept these 
models. 


4.1.3. Porous Grains? 

As a final step in our preliminary study of the dust com- 
position, we consider a grid of models consisting of sili- 
cate grains with different porosities, from compact silicates 
to extremely fluffy silicate grains (V max = 99.375%). The 
marginal probability distribution for the porosity V re- 
veals a bimodal, somehow paradoxical probability distri- 
bution: the most probable models are the ones with either 
no porosity (discussed in Sec. 4,1.1), or very high porosity 
(V = 97.5%). Both solutions are presented in Tab. 6. A 
visual inspection of Figure 8 resolves the ambiguity. The 
compact grain model better reproduces the minimum and 
maximum wavelength ends of the SED, and in particular 
fits the sub-mm point, but it fails at reproducing the inter- 
mediate wavelength range around the emission peak. On 
the other hand, the highly porous grain model is orders 
of magnitude below the observed value at 870 gm, but it 
is within the error bars in the other regions of the SED, 
especially at Herschel wavelengths. In sum, none of these 
models lead to overall good fits to the whole SED. 


4.2. Construction of an Advanced Grain Model 

From what precedes, we conclude that none of the 2- 
cornporient mixtures is sufficient to provide a satisfactory 
fit to the SED. This is essentially a consequence of the disk 
surface density being known a priori in this model; it can- 
not be adjusted to compensate for the emission properties 
of the dust grains. We note, however, that the addition of 
one component to silicates, should it be carbonaceous ma- 
terial. amorphous ice. or porosity, provides some slight im- 
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Porosity 



Fig. 10. Bayesian probabilities for the grain composition in the case of 4-material dust mixtures. 


provements. This suggests that their combined effect could 
produce a more realistic model. 

We evaluate the improvement provided by a 3- 
coinponent mixture on the fit to the SED. We showed in 
See. 4.1.1 that it is not clear how to disentangle the re- 
spective contribution of silicates and organic materials by 
looking at the SED. Instead of including an additional free 
parameter, we consider grains consisting (approximately) 
in the fixed Si/C ratio described in Sec. 3.4 (vc/vsi = 2), 
restricting the discussion to the ice fraction and the porosity 
of the grains. We thus assume that grains are made of sil- 
icate and carbonaceous material, which are either icy ( v ice 
parameter) or porous (V parameter). The results are re- 
ported in Tab. 6. For porous ice-free grains, the best model 
is found for a porosity V = 95% and leads to some im- 
provement to the fit with x 2 — 6.05. If we now consider 
only compact icy grains, we can achieve a. very reasonable 
Xr = 4.48 with comparison to previous models, assuming 
an amorphous ice fraction r ice = 0.7. Overall, both the 
addition of either water ice or vacuum to silicate + car- 
bonaceous grains do contribute to improving the fit to the 
SED of HD 181327. This can can be seen in Figure 8 (bot- 
tom panel). Nevertheless, the fits are still unsatisfactory 
(Xr much larger than 1), and justify an even more advanced 
dust model. 

We thus elaborate a 4-component grain model inspired 
by the “cold coagulation dust model” by Li & Limine 
(2003a, b), which has been successfully employed for mod- 
eling debris disks (Li & Greenberg 1998; Augereau et al. 
1999; Li & Limine 2003a, b; Kohler et al. 2008). We keep 
the Si/C volume ratio fixed (t>s;:t>c=l:2) and let both the 
ice fraction and the porosity vary, in addition to the other 
usual free parameters in this study, namely the minimum 
grain size a m j n , the slope of the size distribution «, and 
the total mass of the disk Af c j ust . This leads to a 5-free- 
parameter disk model to fit the 48 data points in the SED. 
We compute 3 grids of models corresponding to the 3 sur- 
face density profiles shown in Fig. 5 for |<?hgI = 0.30 (i.e. the 
nominal and ±1 a profiles). The Bayesian analysis is then 
applied to the 3 grids separately and the probabilities for 
each point of the parameter space are calculated as the sum 
of the probabilities for the 3 profiles. The detailed results 
are documented in Tab. 6. 

Figure 4 presents the results of the Bayesian analysis on 
the composition of the dust grains, with a 2D-inap show- 
ing the probability as a function of V and t>i ce for the grid 
models, after projection of the probabilities over the other 
parameters, it clearly illustrates the need for a composition 
within a limited range of possible values. In particular, the 


fraction of H 2 0-dommated ices is well defined: the proba- 
bility peak is obtained for an ice fraction V\ ce = 0.67 ±0.07. 
The porosity is not so strongly constrained but we can find 
a peak of probability for V = 0.63 ± 0.21. Eventually, the 
best fit model in terms of x 2 is found for grains made of 
23% carbonaceous materials, 12% silicates, 65% ice, and 
with a porosity V = 65%. As can be seen in the bottom 
panel of Figure 8, for such grains, the synthetic SED very 
reliably describes the observed SED from the mid-IR to 
the sub-nun, with yf = 1.64, a value significantly better 
than those obtained for simpler grain models. The grain 
size distribution that comes out. of the statistical analysis 
has a m ,„ = 0.81 ± 0.31 fan and k = —3.41 ± 0.09, a power 
law r index very consistent with that predicted for theoretical 
collisional cascades (k — —3.5). The mass in grains smaller 
than 1 mm, Md ust — 0.051 ± 0.016 M®, may look surpris- 
ingly small at first glance compared to the other models, 
but it results from the significant fractions of ice and vac- 
uum in the grains, implying larger mass opacities (Fig. 3.4). 
The temperature of the grains at the surface density peak 
position (89.5 AU) ranges from 86 K for the smallest grains 
to 40 K for ~ 8 mm-sized grains. 

To summarize, our approach, which is based on a fixed 
surface density profile constrained by resolved imagery of 
the disk, clearly highlights the need to consider the presence 
of both ice and porosity to find an overall good fit to the 
SED. Icy grains are naturally expected at such a large dis- 
tance from a solar-type star and the fluffiness of the grains 
is likely a by-product of their formation process. We ap- 
plied a similar cold-coagulation dust model as Li & Lunine 
(2003a) did to fit the SED of the very young HD 14 1569 A 
debris disk, except that we used no a priori assumption 
on the ice fraction. By arguing the grain composition must 
reflect the primordial abundances of condensable elements, 
they chose v ice « 30% for cold regions and they find a best- 
fitting porosity V ss 0.8. Hence we propose a dust model 
that is especially ice-rich with a degree of porosity that, is 
not unexpected. Further comparisons with the literature 
and implications on the nature of the dust are discussed in 
more details in Sec. 6. 


5. Limits on the Gas Content 

Aged 12 Myr-old, HD 181327 offers a chance to assess the 
amount of gas in a young debris disk, at a time when, ac- 
cording to planet formation theories, the accretion of gas 
by giant planets has already occurred while the formation 
of rocky planets might be ongoing. 
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We attempt to use the upper limits on the oxygen and 
ionized carbon fine-structure lines ([O I] at 63 and 145 pm, 
and [C II] at 158 //in. see Tab. 3) to constrain the amount of 
gas present in the disk. We solve the photochemistry and 
energy balance in the disk using the photochemical code 
ProDiMo 2 . The abundance of 71 species is computed in 
steady-state, under the influence of adsorption, desorption 
(thermal, photo and cosmic-ray induced), and molecular 
hydrogen formation on grain surfaces. 

The gas density profile is supposed to match the one 
of the dust, as presented in Sec. 3.1. The total mass of 
solid (defined as the mass in grains from 0.81 //m to 8 mm) 
is constrained by the best dust model of Sec. 4.2 (M so iid 
= 4.9 x 10~ 7 M© = 0.164 M$) and we let the total 
atomic+molecular gas mass A/ gas vary from a gas-poor disk 
( AT gas / M S oiid = 1 0 ~ 2 ) to a gas-rich disk (Af gas /M so i k i=10 :i ). 
The parameters of the model are summarized in Tab. 7. 
One of the main parameters controlling the gas energy bal- 
ance is the amount of PAHs (/pah), which remains uncon- 
strained in the disk of HD 181327 because tlieir abundance 
is too small to produce detectable PAH emission features in 
the IR. We adopted two ’’extreme” values for /pah: a low 
PAH abundance (/pah=10~ 5 * ) and a high PAH abundance 
(/pah=0.1 * ) 3 . 

The non-LTE line fluxes for the [O I], [C II], CO J=3 
-* 2, and CO J— 2 -¥ 1 transitions are given in Fig. 11 and 
Appendix B. Among all the lines, only the CO J= 3 -> 2 
and J= 2 — > 1 line fluxes behave monotonically, increasing 
in flux with increasing disk mass. The reasons for the appar- 
ent erratic behaviour of the fine-structure lines come from 
the non-linear nature of several physico-chemical phenom- 
ena. For example, the [C II] flux depends on the abundance 
of C + and on the excitation of the first energy level. The C + 
abundance first increases with total gas mass until CO self- 
shielding against photodissociation is efficient enough such 
that a significant fraction of C + is converted into neutral C 
and CO. The presence of CO also changes the heating and 
cooling balance, since the CO absorption of IR stellar pho- 
tons becomes the main source of heating. The change from 
a disk with no CO self-shielding to a CO self-shielded disk 
manifests itself as a jump in CO fluxes between a disk with 
gas-to-solid ratio of 1 to a disk with gas-to-solid ratio of 10. 
The abundance of PAHs has a weak influence on the gas 
temperature because the small UV flux from a F5/F6 star 
without excess accretion luminosity cannot ionize PAHs. 
On the other extreme, when the gas mass is very small 
(gas-to-solid mass ratio below 1), the PAH and molecular 
absorption heating are nonexistent. Photoelectric heating 
by large grains dominates and explains why the line fluxes 
are similar for the high and low PAH abundance models. 

Line fluxes from models with /pah = 10 -5 are consistent 
with the PACS upper limits for all gas-to-solid mass rat ios 
(upper part of Table B.l). Only in the case /pah=0.1. the 
disk gas mass is constrained by [01] at 63 pm to be below 
~17 M© (iWg as /.M so iid < 10 2 , lower part, of Tab. B.l). 

Eventually, even if the integration of the [O I] 63 and 
145 pin, and [C II] lines we performed was not sufficient 


2 Woitke et al. (2009) provide a detailed description of the 

code, with subsequent improvements described in Kamp et al. 

(2010), Thi et al. (2011), and Woitke et al. (2011). 

J The high PAH abundance is an upper-limit, to the amount 

of PAHs in HerbigAe disks, where it, is found to be of the order 

of 0.1-0.01 (e.g., Meeus et al. 2010, /pah =0.03 for HD 169142). 


to derive strong constraints on the gas-to-solid mass ratio 
in the disk, the combination of the GRaTeR and ProDiMo 
codes proved efficient to make predictions on the fluxes in 
the CO /= 3 — i 2 and J — 2 -A 1 lines. These lines will be 
detectable for instance by ALMA, finally allowing to assess 
the gas content of the debris disk. 


Table 7. Setup of ProDiMo for the gas disk model 


Parameter 

Value 

Elemental abundances 

ISM-like 

UV field 

1 Habing field (a) 

Cosmic-ray ionization rate Q 

1.7 x 1(U 17 s” 1 (b) 

Non-thermal line width dv 

0.15 kin.s -1 

/pah ! c) 

10^ 5 - 10’ 1 


^ Dalgarno (2006), (l,) Standard interstellar UV field. (c) /pah 
scales with the gas mass and is defined with respect to an inter- 
stellar abundance of 3 x 10~‘. 



0.01 0.10 1.00 10.00 100.00 1000.00 
Gas-to-solid mass ratio 


Fig. 11. Expected gas line fluxes as a function of gas-to-solid 
mass ratio for either a low-PAH abundance (thin lines) or a 
high- PAH abundance (thick lines). The horizontal lines denote 
the PACS 3tr upper-limits and do not include the 30% flux cali- 
bration uncertainty. The gas-to-solid mass ratios assume a solid 
mass of 0.164 M®. 


6. Additional constraints on the disk 

6.1. Scattered light: evidence for irregular aggregates 

The surface density profile used in this study comes from 
a new reduction of the scattered light images of Schneider 
et al. (2006) obtained with HST/NICMOS (see Sec. 3.1). 
This ensures by definition that the shape of any synthetic 
scattered light surface brightness profile in this study per- 
fectly agrees with the observed profile at 1.1 pm. However, 
the consistency is limited to the shape, and the agreement 
in terms of absolute flux is not guaranteed. This is because 
our fitting approach in Sec. 4 consists of adjusting the disk 
mass to get a best, fit to the SED, which is equivalent to 
scaling up or down the surface density profile. A reason 
for adopting this strategy comes from the fact that the 
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scattered light images do not directly provide the surface 
density E(r), but its product by the mean scattering cross 
section <7 sca at the observing wavelength (1.1 fim). The fit- 
ting process to the SED, on the other hand, yields an in- 
dependent,, absolute scaling of the surface density profile, 
thus allowing one to obtain a first estimate of <r sca . We find 
cr sca = 3.7/un 2 in the NICMOS 1.1 gm band. 

A second estimate of <r sca is obtained by calcu- 
lating the theoretical scattering cross section for com- 
pact spherical grains given the grain properties (com- 
position, size distribution) that were inferred in Sec. 4: 
<r S ca(A) *= 7ra 2 <3sca(a, A)dn(a) where Q SC[l is the 

scattering efficiency at wavelength A and for a grain radius 
a. For the best model (in the Bayesian sense, Tab. 6) we 
find <r sca = 16.8 gm 2 at 1.1 //m, a factor of about 4.5 higher 
than the first estimate. This implies a peak flux density of 
3.14 mJy.arcsec -1 inconsistent with the 0.76 m.Jy.arcsec™ 1 
measured in the NICMOS profile. In other words, this 
means that our model predicts an albedo (ui = 59% at 
1.1 /an) that is a factor of about 4.5 larger than the ob- 
served value. This suggests that smaller grains than those 
required to fit the SED dominate the scattering behavior 
of the disk. It is worth noting that such a conclusion is not, 
specific to HD 181327, but is also observed in other debris 
disks. Recent examples include HD 207129 for which Krist 
et al. (2010) predict an albedo (based on dust properties 
also obtained from a fit to the SED) that is an order of 
magnitude larger than the value suggested by the scattered 
light image. A similar inconsistency is noticed for HD 92945 
by Golimowski et al. (2011) who point out that different ob- 
serving diagnostics suggest different minimum grain size. 

It is important to keep in mind that although our model 
suggests the grains are very porous and hence likely to have 
an irregular, fractal structure, we used hard spheres (Mie 
theory) to calculate both the thermal emission and scat- 
tered light optical properties. Unlike the thermal behavior, 
light scattering may strongly depend on the exact shape 
of the dust aggregates, and on the size of their elemental 
constituents, two features that are not accounted for by 
the Mie theory (Voshchinnikov et al. 2007). Another piece 
of evidence that Mie theory might be a limiting factor for 
a self-consistent model of HD 181327 comes from the ob- 
served scattering asymmetry parameter {/hg that is a mea- 
sure of the direction of the light scattered by a particle. 
Our best model (Sec. 4.2) has a mean scattering asymme- 
try parameter |<?hg|= 0.95 at A = 1.1 /an, much larger than 
the \(/hcI— 0.30 measured on the NICMOS image. Only a 
population of much smaller particles could produce such a 
low scattering asymmetry parameter if we stick to t he Mie 
formalism, as already noted by Schneider et al. (2006). An 
analogous problem arises from the modeling of the debris 
disk of HD 92945 (Golimowski et al. 2011). Our micron- 
sized grains are highly anisotropic, essentially due to their 
spherical shapes. Incorporating vacuum homogeneously in- 
side the grains does not suffice to decrease [<?hg| and a 
more accurate description of the shape of fluffy aggregates 
is eventually needed. In fact, the scattering efficiencies and 
asymmetry parameters of fluffy aggregates are known to be 
correctly described by an EMT-Mie model (porous spheres) 
only in the Rayleigh regime, i.e. when the size of the in- 
clusion (vacuum here) is small compared to the wavelength 
(see e.g. Voshchinnikov et al. 2007). This is not the case 


anymore for the relatively large grain sizes and porosity at 
stake when the wavelength becomes small. 

6.2. A narrower parent belt? 

Unlike scattering, we expect the continuum thermal emis- 
sion to be largely independent of the exact shape of the 
grains. Our best model of HD 181327’s disk produces a 
sub-mm flux of 39mJy at 870 gm, consistent within 2o 
with the LABOCA measurement (52±6mJy). The ATCA 
measurement at 3.2 mm was not model-fitted, because of 
large uncertainties on the optical constants at millimetre 
wavelengths. As a sanity check, we calculate the flux in the 
ATCA band using a power-law extrapolation of the SED. 
In the far-IR, to mm regime, the model emission falls off 
as ^-2-51*0.164 xiiis results in a flux equal to ITIq 4 niJy, 
which is within ~ 3cr compatible with the measured flux 
of 0.72 ± 0.25 in Jy. Thus, the model remains valid for the 
longest wavelength observations (coldest emitting parti- 
cles), showing that the ring resolved in the near-IR with 
the HST is compatible with the spatial distribution of the 
largest grains, at least at first order. 

A limit of our approach is the assumption that the grain 
size distribution is the same everywhere in the disk, while 
it is known that only the smallest grains can populate the 
outer disk due to radiation pressure. The maximum size 
and power law slope of the dust size distribution should 
vary with the distance from the star (e.g. Fig. 4 in Augereau 
et al. 2001), whereas the model keeps both large grains in 
the outer disk and a constant slope for the size distribution. 
Such a radial variation of the size distribution is suggested 
by a comparison of a synthetic image at 70/mi with the 
PACS 70/nn map (Sec. 2.1). The synthetic image appears 
broader than the observed one with a FWHM of 7.08". 
versus 6.44" for the PACS image. Our interpretation is that 
the PACS band reveals larger grains that experience less 
the force due to the radiation pressure and thus are more 
bound to the parent belt. Assuming all grains are equally 
distributed results in putting too many large grains in the 
outer disk: a more complex model should account for the 
true dynamics of the particles. 

6.3. Emission from inside the belt? 

The region of the disk interior to 75 ALT cannot be reliably 
assessed from the NICMOS data because of the corona- 
graph and artefacts from the PSF-substraction below 1.5" 
(78 AU). Therefore we cannot exclude the existence of an 
additional dust component interior to the main belt, that 
would be invisible in the surface density profile we use. If 
micron-sized grains populate the inner disk, they would ra- 
diate mainly in the mid-IR. which might affect the results of 
our study. From the Gemini/T-ReCS images of Chen et al. 
(2008). we estimate that less than ~15% of the light emit- 
ted by the disk at 18.3 /mi originates from a region interior 
to 1" (Fig. 2 in Chen et al. 2008). 

We perform a simple test of the robustness of our re- 
sults. Suppose that a belt of black body dust particles is 
orbiting the star at a distance r < 1". We calculate a black 
body spectrum than peaks at 18.3 pm (T = 150 K) and 
we adjust it, such that, its maximum flux represents 1 5% of 

4 Note that, at longer wavelengths, the emission seems to falls- 
off faster than a simple power-law. 
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the 53mJv attributed to the disk at, that wavelength. Then 
we subtract the black body flux from the observed SED 
of Tab. 4 to simulate an observation of the system exclud- 
ing the hypothetical inner component. We reproduce the 
approach of Sec. 4.2 with this new SED. Interestingly the 
results of the Bayesian analysis are barely impacted. To ac- 
count for the deficit of flux at mid-IR wavelengths, only the 
minimum grain size needs to be slightly adjusted to elim- 
inate the smallest grains: a m i n = 1.0 pm and k = —3.43, 
while we have a m ; n = 0.89 pm and n = —3.42 for our best fit 
(Tab. 6). The composition remains unchanged for the best 
model (x, 2 = 1.52) and no change is induced that would im- 
prove the consistency of the model with the scattered light 
flux (Sec. 6.1). Consequently we cannot reject the possible 
existence of an inner component to the disk, but none of 
our conclusions would suffer from its existence. 


7. Grain dynamics and ice origin 

7. 1 . Blowout size 

Grains in the HD 181327 debris disk for which radiation 
pressure overcomes gravity are put on unbound orbits and 
eventually ejected from the system on a short time scale 
(see for instance Krivov 2010, and ref. therein). The balance 
between radiation pressure (Fp) and stellar gravity (Eg) 
forces is evaluated with the distance independent, but. .size- 
dependent, 8 pr factor, given by 

o _ ]£H _ a (Qpr) L * (2 \ 

Ppl ~ \F a \~ mATrcGM* 

where (Q pr ) is the pressure efficiency averaged over the 


Si 

Si+am.ice 



Fig. 12. 0 pT = |Ep|/|Fg| as a function of grain radius 


stellar spectrum, m the grain mass, a the grain geometri- 
cal cross section, c the speed of light, G the gravitational 
constant, and £* and Af* documented in Tab. 1. Assuming 
that the grains are spherical, Eq. 2 yields a critical size be- 
low which the grains originating from the parent belt are 
ejected from the system: flbiow — a(0 p r = 1/2) (assuming 
circular orbits for the parent bodies). 

In Fig. 12. we represent the /3 pr ratios as a function 
of grain radius for different compositions. Our best 4- 
component grain model has a blowout size flbiow = 4-83 pm, 


that is 5 times larger than the inferred minimum grain size 
in the system. Models implying no porosity nor carbon, 
on the other hand, have blowout sizes of the order of the 
inferred value of a min , more consistent with the theoreti- 
cal argument that grains below the blowout size should be 
expelled from the system. The response of silicate grains 
to the stellar irradiation and gravity fields is only slightly 
affected by the presence of ice, leading to similar shape 
and moderated values for /3 pr (a). Organic refractory grains 
on the contrary are highly sensitive to the radiation pres- 
sure and they can reach high /3 pr values. Fig. 12 shows that 
porosity has the effect of flattening the /3 pr («) function, 
pushing the critical value /3 pr =1/2 to larger grains. Hence 
the addition of porosity to a mixture of carbon and sili- 
cates produces a strong effect: ;3 pr not only reaches a high 
value at its maximum due to carbon, but the /3 pr (a) func- 
tion decreases slowly such that the critical value /3 pr = 1/2 
is reached only for large grain sizes. 

Decreasing the amount of carbon in the grains which 
was imposed somehow arbitrarily in our model (see Sec. 4.2) 
would mechanically decrease their blowout size. Pure sili- 
cate grains with the same porosity and ice content have 
flbiow = 0.89 pm, but they do not provide good fits to the 
SED. We computed additional grid of models with alter- 
native vc/vsi values. The results confirm the conclusion of 
Sec. 4.1.1: all Si+C-fHaO+porosity models lead to similarly 
good fits to the SED, as long as a sufficient fraction of car- 
bon is considered. For instance, if silicates represent only 
10% of the C+Si mixture, then v- lce remains at 0.65 and 
the porosity is reduced to 55%, with A/dust = 0.043 Af®, 
a m i n = 0.89 pm and \f = 1.66. All these values remain 
within the error bars of the original model. Models involv- 
ing carbon-poor grains do not provide good fits to the SED, 
but as long as «cA'Si > 0.5, the results remain unchanged. 
In summary, the discrepancy between a mm and <ibiow does 
not lie on the hypothesis we made on the composition of 
the grains. 

Our estimation of Ubiow is based on the assumption that 
the grains are spherical and homogeneous, but the actual 
response of inhomogeneous aggregates to radiation pressure 
cannot be that simple. Saija et al. (2003) for instance show 
that the fluffiness of the grains tends to decrease their ra- 
diation pressure cross-section, as if aggregates behaved like 
their elemental constituents. This would make grains below 
the theoretical blowout size less sensitive to radiation pres- 
sure, and allow them to remain bound to the star. Thus the 
apparent inconsistency can be seen as another signpost of 
the complex structure of the grains, and we highlight that 
fixing a„,i n as equal to «biow calculated with the (EMT-)Mie 
theory would not be a relevant approach when dealing with 
porous grains. 

Another possibility is tiiat, we are observing a popula- 
tion of grains that is not in equilibrium (i.e. an imbalance 
between production and ejection), following for instance a 
major collisional event, such as an asteroid break-up. An ex- 
cess of dust grains from all sizes would have been produced, 
including grains below the blowout size. These grains could 
have survived destructive collisions and radiation pressure 
blowing for long enough to be observed. Of course, observ- 
ing such an event is highly unlikely, because dust typical 
lifetime is only of the order of 1,000 years, which is small 
with comparison to the age of the star (see Sec. 7.2). But 
this “rare event hypothesis” would corroborate the scenario 
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proposed by Chen et al. (2008) to justify the brightness 
asymmetry they claim to observe in thermal light. 

7.2. Time scales and origin of the ice 

Collisions occurring between grains in a debris disk are 
a key process as it generally results in the production of 
smaller dust particles, either by fragmentation of the tar- 
get, or by its excavation. Any phenomenon having a longer 
characteristic time scale than the collisional time scale will 
be considered negligible in the evolution of the dust popula- 
tion. Based on the principle that any particle will cross the 
disk midplane twice during one orbit, Thebault & Augereau 
(2007) proposed a simplified relation for the collisional time 
scale, valid for the smallest grains: t m \ — torbitAx 80 ) w hh 
f-orbit the orbital period at the surface density peak, and rJ eo 
the geometrical vertical optical depth. Here the orbital pe- 
riod at 89.5 AU is ~ 850 years, leading to f co j ~ 7600 years 
for the best model. This is more than 3 orders of magni- 
tude smaller than the age of the system and 2 orders of 
magnitude smaller than the characteristic time-scale for 
Poynting- Robertson drag: t pr = 0.1 Myr for j3 pr = 1/2 
grains (Wyatt 2008). The evolution of the HD 181327 de- 
bris disk is clearly dominated by collisions and it requires 
a mechanism to generate dust on a 10 3 — 10 4 years charac- 
teristic time scale. 

Attempts to model time-dependent chemistry with 
ProDiMo reveals that the grain water ice mantle shall be 
photo-evaporated in a few 100 years (for ~ 1 pm grains), to 
a few 1000 years (for ~ 10 pm grains), a similar (or shorter) 
time scale to the dust mutual collision, requiring as well an 
ice reprocessing mechanism. This provides a consistent view 
that the icy dust grains are maintained through collisions 
in a reservoir of icy planetesimals. We highlight that, while 
detecting the mid-infrared spectral features in the emission 
of a debris disk is very challenging - unambiguous detec- 
tion have only been made in brighter circumstellar disks 
(see, e.g. Malfait et al. 1999) the analysis of the contin- 
uum SED alone can provide clues on the amount of ice in 
dust disks, thus on the nature of the parent bodies, often 
referred to as the leftovers of planet formation. 

In the Solar System Kuiper Belt, the study of icy plan- 
etesimals, Kuiper Belt Objects (IvBOs) in particular, is 
a prolific topic (for a recent review see e.g. Jewitt 2010). 
The presence of ice, together with the porous nature of 
the objects, is inferred from estimations of their densities. 
Unambiguous detection of water ice has been made pos- 
sible by the spectroscopic near-IR signature of crystalline 
ice, for instance by Delsanti et al. (2010) who complete 
their study of the plutino Greus by radiative transfer con- 
siderations to conclude that it contains a mixture domi- 
nated by amorphous and crystalline water ice. Confronted 
by space weathering, crystalline ice should be amorphized 
in a short timescale. This implies that it needs to be stored 
(and possibly produced) in the interior of the object before 
being brought to the surface, possibly through eryovolcanic 
events. The nature of water ice in KBOs remains unclear, 
but the vast majority of these objects are too small to ex- 
perience such geological processes (Dumas et al. 2011, and 
references therein). Thus KBOs covered by crystalline ice 
can only constitute particular examples, and the water ice 
covering the collisionally produced dust in the Solar System 
is more likely amorphous. The youri'g age of IID 181327 al- 
lows the possibility that primordial crystalline ice has sur- 


vived in its debris disk. Our study favors the amorphous 
ice, suggesting that either ice was predominantly incorpo- 
rated in an amorphous form in HD 181327’s Kuiper Belt 
objects, or that the amorphization timescale is shorter than 
~10 Myr in this distant belt. 

8. Discussion 

The properties of the HD 181327 debris disk can be com- 
pared to those of the disks surrounding two of its co- 
inoving ’’sister” stars, 8 Pictoris (A3V, 1.7 Mq) and 
AU Microscopii (M1V, 0.5 Mq), which formed in the same 
molecular cloud from the same material. Of course compar- 
isons are limited by the specific models used to study each 
object. The relatively large porosity found here should be 
put in the context of the grain properties found for AU Mic 
(Graham et al. 2007: Fitzgerald et al. 2007). Assuming a 
Si+C+H>0 mixture, only highly porous grains (V = 80%), 
allow to reproduce altogether the scattered light profiles 
(exhibiting strong forward scattering), the high polariza- 
tion and the thermal SED. In that model, the Mie theory 
remains efficient to represent the scattering properties of 
porous grains, because the grains filling the outer disk are 
very small and stay in the Rayleigh regime. A total dust 
mass of 0.01M® is found, that essentially reflects the con- 
tribution from the inner mm-sized grains, concentrated in 
a ~40 AU belt. A similar porous grain model was used by 
Li & Greenberg (1998) to model 8 Pictoris yielding a dust 
mass ~0.4 M®, contained in the ~ 100 AU-wide outer disk. 
Overall, ’’cold coagulation-like” dust models prove effective 
at reproducing observations of these young debris disks. A 
trend appears for increasing dust mass with spectral type 
as could be intuitively expected, but despite this difference 
of scale, the three stars are orbited by a similar material. 

Henceforth, given the young age of these objects, we 
are here investigating the conditions holding in disks just 
after the planetary embryos may have formed, the gas dis- 
persed. and possible accretion onto outer giants happened. 
To explain the belt-like structure of the HD 181327 disk, 
the most straightforward explanation is to invoke the dy- 
namical influence of a planet inside the belt. The Wisdom 
( 1980) analytical criteria gives the width of the chaotic zone 
a planet mean-motion resonances would create. Assuming 
circular orbits for the planetesimals, as suggested by the 
sharpness of the belt, this yields a range of possible val- 
ues for the planet semi-major axis depending on its mass. 
Assuming the inner edge of the disk is located at 78 AU 
from the star, the yet undiscovered planet would lie from 
73 AU for a Neptune-mass planet, down to 59 AU for a 10 
Jupiter-mass planet . 

The material at that stage appears as very hydrated 
providing a considerable reservoir of water in the outer 
system. The existence of such a reservoir is hypothesized 
in scenarios of the Solar System formation to explain the 
unexpectedly large water content on Earth for instance 
(Morbidelli et al. 2000). Indeed, terrestrial planets at this 
age (12 Myr) are likely too hot to hold on to water and 
later impacts of icy comet-like bodies from the outer re- 
gions are a possible mechanism for delivering water and 
volatiles to planets in the habitable zone. Thus, not only 
HD 181327’s debris disk reveals that the system has already 
grown planetary embryos including maybe the bricks of 
rocky planets in the inner disk but it likely holds the con- 
ditions that were prevailing before the Solar System Late 
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Heavy Bombardment that depleted the Kuiper Belt and 
determined terrestrial planets history. 

9. Summary and conclusions 

We present new Hersehel/PACS and ATCA imaging of the 
HD 181327 debris disk, a 12Myr old star belonging to the 
p Pictoris moving group. The 70 and 100 //m images as 
well as the 3.2 millimetre observation reveal a marginally 
resolved extended emission. Fluxes in the three PACS pho- 
tometric channels were derived and complemented by con- 
tinuum fluxes obtained from PACS spectroscopic obser- 
vations. They provide precise information on the contin- 
uum thermal emission of the debris disk, but detailed mod- 
els such as the one developped in this study are required 
to interpret the precise nature of the emitting material. 
Degeneracies between the grain properties and disk’s spa- 
tial structure are broken in our model by making use of 
a reprocessed HST/NICMOS coronagraphic image reveal- 
ing that the dust is bound to a narrow belt peaking at 
89.5 AU from the. central star. We thus focus our study 
on the dust properties (composition, size distribution) us- 
ing a radiative transfer code, and run a large grid of de- 
bris disk models for subsequent statistical (Bayesian) anal- 
ysis. We conclude that the thermal light, spectrum of the 
disk is well reproduced assuming a mixture of silicates, 
carbonaceous material and amorphous ice, as well as vac- 
uum to mimic a significant, porosity. The observed SED 
could be very reliably reproduced by a population of grains 
larger than ~ 1 /mi, which follows a power-law distribution 
matching the ideal collisional equilibrium case and repre- 
sents ~0,05AT© (1.5 x 10 7 A / . } in grains with sizes up to 
1 mm. Extending the size distribution to the parent bodies 
of the observed dust grains leads to a mass of ~ 50 AT© for 
planetesimals up to 1 km. The composition we infer for the 
dust grains must be a good representation of a considerable 
reservoir of material hidden in unseen planetesimals. 

In this study, we attempted to use Herschel/PACS spec- 
troscopy to detect for the first time far-infrared emission 
lines in the debris disk. The non-detections we obtain were 
used to investigate the gas content of the debris disk with 
our photochemical disk model ProDiMo based on the dust 
model found with GRaTer. The main uncertainty lies in the 
amount of PAHs in the disk. The current non-detections 
of [OI] and [CII] lines alone do not provide unambigu- 
ous upper limits to the gas content. Only in the case of 
a high PAH abundance, we can set an upper limit, to the 
gas mass of ~ 17 AT© (one Neptune mass). Coupling with 
other tracers, in particular with the CO lines accessible to 
ALMA, offers much better prospects to reach lower limits 
on the low gas content of debris disks. However, steady- 
state chemistry is likely not ideally suited to model debris 
disk, and in future studies more appropriate prescriptions, 
time-dependent chemistry in particular, will be explored. 

We note that both the low albedo estimated from the 
scattered light images and the dust dynamics are not well 
modelled by spherical grains filled with vacuum, and that 
the grains are likely complex fluffy aggregates. Though their 
thermal behavior is correctly described by the EMT-Mie 
theory, the scattering of light at short wavelengths and 
their sensitivity to radiation pressure require a more accu- 
rate description of the exact shape of the aggregates. Our 
statistical analysis highlights that ignoring the porosity of 
the grains or setting the minimum grain size as equal to 


the blowout size would result in incorrect interpretations of 
the observations, in particular because the grain opacity is 
strongly dependent on the exact grain composition. 

We stress some necessary limitations of the models. 
First, we use a grain size distribution that is independent 
of the distance to the star, whereas it is known that only 
the smallest grains can populate the outer parts of the disk 
under the effect of radiation pressure. In this paper we dis- 
cussed clues that the observed scattered light profile may 
not be a perfect representation of that seen in thermal light, 
by Hersehei or ATCA. The long wavelength emission em- 
anating from the largest grains might originate from an 
even narrower belt. Further studies will need to account, for 
the radial differentiation of grain sizes with the objective 
to fit altogether resolved images from the near-IR to the 
millimetre domain. Secondly, the Dohnanyi-like power-law 
size distribution is known to be an insufficient description; 
if tiie real distribution implies wavy patterns as proposed 
by Thebault & Augereau (2007), then adjusting a power- 
law likely puts too many/too few small grains in the disk. 
Dynamical models are required to understand better the 
general outcome of the collisional evolution of the disk. 
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Appendix A: Statistical analysis of the results 

The strength of the GRaTer code, in comparison with an 
optimization algorithm that would provide a unique best- 
fit, solution and possibly miss local minima, is that all the 
models are stored enabling fast and easy post-processing of 
the results. 

To find the best solution in the grid of models, we use 
a classical least-square minimization approach with ttf ree = 
3 to 4 free parameters depending on the grain composition 
and N = 48, which results in 43 to 45 degrees of free- 
dom (dof) (Tab. 6). The free parameters are k, a m in and 
-Mdust, and depending on the situation «c, v- lce or V . Mdust 
is handled separately because it is automatically adjusted 
to reproduce the surface density profile in the code. 

Let a be the vector of parameters. Here 

a = (k, a m in, f ice, L 5 ). For each set, of parameters, 

we can define the quantity: 


E 


(A - Mi ) 2 


(A.l) 


where A and AT,- are the observed and modeled values, 
respectively, and a,- is the measured uncertainty. For each 
X 2 value, we can define the likelihood of the data, i.e. the 
probability of the data given the parameters of a model. 
This assumes the errors are distributed normally and are 
uncorrelated. 

P(D\a) = ~exp(^^j (A. 2) 
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where C is a normalization constant and D = {jD('i)}i=i...N. 
In order to evaluate the effect of each individual parame- 
ter tti on the quality of tlie fit, we implement a Bayesian 
inference method. We define the posterior probability of a 
as the probability distribution of the parameters a given 
the data D. It has to be noted that, this formulation im- 
plicitly assumes that the model is suited to describe the 
observation. Bayes’ theorem says that this quantity can be 
written: 


P(a\D) 


P(a)P{D\a) 

P(D) 


(A-3) 


where P(D) is the prior probability of the data under all 
possible hypothesis (a normalization constant in practice), 
and P(a) is the prior probability distribution of the param- 
eters a. Here we use uniform priors for all the parameters, 
meaning that we do not favor any realization of the parame- 
ters prior to obtaining the data. That fundamental theorem 
provides a simple way to independently assess the param- 
eters of a model. Indeed, the marginal probability of each 
parameter a,\ can be computed regardless of the values of 
the other parameters a^\: 


P{a\\D) = J... J J ... j P(a^\D)da 0 ...da ]tm 

«o a*-! a »+ 1 a, i m ax 

(A.4) 

Then we obtain a probability law as a function of the re- 
alization of each individual parameter. Finally, we estimate 
the parameters by calculating the mathematical expecta- 
tion and variance. These are the values presented in Tab. 6 
where the uncertainties are 3rr. 




E(a,i) = J aidP(a-,\D) 

(A.5) 

a 'hnin 

<y 2 (ai) = E(af) - E(a ;) 2 

(A.G) 
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Table B.l. Gas line fluxes. The gas to solid mass ratios assume a solid mass of 0.164 \L. . The PACS 3-cr upper limits do not 
include the 30% flux calibration uncertainty. 


Gas-to-Solid 

[O I] 63/nn 

[O I] 145/rm 

[C II] 158/im 

CO .7=3 -a 2 

CO ,7=2 -> 1 

mass ratio 

(W m“ 2 ) 

(W nr 2 ) 

(W m“ 2 ) 

(W m“ 2 ) 

(VV m~ 2 ) 


Herschel-PACS (3a) 




9.8xl0“ 18 

8.5 xlO -18 

7.0xl0“ 18 

- 

- 



ProDiMo predictions with /pah=10 


1000 

2.9xl0“ 18 

1.6x10“ 18 

2.6xlO“ 20 

1.7xl0“ 19 

•5.6xlO“ 20 

100 

6.7x10“ 19 

1.2xlO“ 20 

4. lx 10“ 20 

1.5xl0“ 19 

5.3xlO~ 20 

10 

2.1xl0“ 18 

8.2x 10“ 21 

2.2xl0“ 19 

1.3xl0“ 19 

3.8 x 10~ 2 ° 

1 

2.4xl0“ 18 

4.0x 10“ 2 ° 

6.2xl0“ 19 

3.4xl0“ 22 

6.1 xlO -23 

0.1 

6.4xl0“ 19 

5. Ox 10“ 2 ° 

1.9xl0“ 19 

l.OxlO -23 

1.6x 10“ 24 

0.01 

2.8xlO“ 20 

2.9xl0“ 21 

1.5xlO“ 20 

1.1 x 10 “ 24 

2.7xl0“ 25 



ProDiMo predictions witli /pah =0.1 


1000 

2.7xl0“ 17 

6.4x10“ 18 

4.1xlO“ 20 

4.8xl0“ 19 

1.6xl0“ 19 

100 

7.1xl0“ 18 

2.0xl0“ 19 

6.0xl0“ 20 

3.1 xlO” 19 

l.Ox 10“ 19 

10 

3.2xl0“ 18 

1.5xlO“ 20 

1.4xl0“ 19 

1.4xl0“ 19 

3.9xlO“ 20 

1 

2.4xl0“ 18 

4.0x10“ 2 ° 

4.2xl0“ 19 

2.6xl0“ 22 

4.6xl0“ 23 

0.1 

6.3xl0“ 19 

4.8x l()“ 20 

1.8x10“ 19 

l.OxlO -23 

1.5xl0“ 24 

0.01 

2.8x 10” 2 ° 

2.9x 10“ 21 

1.5xlO“ 20 

l.lxl0“ 24 

2.7xl0" 25 
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